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NCC 2-585 


ANNUAL REPORT 
Technical Proposal 

The objective of this work is to develop, verify and incorporate the 
baseline two-equation turbulence models which account for the effects 
of compressibility into the three-dimensional Reynolds averaged Navier- 
Stokes (RANS) code and to provide documented descriptions of the mod- 
els and their numerical procedures so that they can be implemented into 
3-D CFD codes for engineering applications. 

Computational Fluid Dynamics (CFD) is recognized as a significant 
engineering design tool in modern engineering projects. The design of 
airplanes and aerospace vehicles is a highly complex process. One of 
the critical tasks involved in such a process is the ability of the turbu- 
lence model to predict boundary layer separation, shock- wave/turbulent 
boundary layer interactions, viscous/inviscid interactions, transition, ad- 
verse pressure gradient flows, rotation and streamline curvature effects, 
mixing problems, and a range of other turbulence phenomena. 

Turbulence models are developed on the basis of insight gained from 
experimental and theoretical research. The complexity of turbulence re- 
quires that the mathematical models be guided by the flow physics in a 
rational and practical approach. Test and validation of new models with 
data of recognized quality is an essential step toward model acceptabil- 
ity. A database with a detailed experimental information on turbulent 
flows is necessary as a benchmark for testing and validation of model 
modifications. The simulations of complex turbulent flows with the new 
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turbulence models complement and assist the design of the data base for 
turbulence modeling validation. 

A major consideration throughout the research effort is the develop- 
ment of improved corrections to the turbulence models and the identifi- 
cation of models which are superior in their predictive capabilities. 
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Work Statement 


The work statement for this period is: 

1. Prepare material and documentation for AIAA paper 94-1905 to be 
presented in the 1994 Applied Aerodynamics conference, to be held 
June 20-23 in Colorado Springs, Colorado. 

2. Supervise and deliver information needed to implement compressible 
two-equation turbulence models into 3-D codes. Prepare material 
and documentation for AIAA paper 94-2950 to be presented in the 
30th AIAA/ASME/SAE/ASEE Joint Propulsion Conference to be 
held June 27-29 in Indianapolis, Indiana. 

3. Study and continue analysis of 3-D turbulence model corrections 
with numerical simulations of 3-D intersecting shock- waves/ turbu- 
lent boundary layer interaction (SWBLI) flows at Mach 8.3 and 4. 

4. Prepare material and documentation of the numerical simulations 
of the 3-D Reynolds Averaged Navier-Stokes equations of the 3-D 
Intersecting Shock Wave - Turbulent Boundary layer Interactions 
(ISWBLI) at Mach 8.3 and 4. 
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Work Accomplished 

The work accomplished is presented with a general description of 
the present general capabilities and the achievements in this performance 
period. These achievements are published in AIAA papers whose general 
descriptions are shown below. 

1. Present General Capabilities. 

The RANS 3-D Reynolds Averaged Navier-Stokes code has presently 
the capability to predict turbulent flows with the following models: 

• Wilcox k-u model (baseline model) 

• Menter BSL k — u / k — e model 

• Menter SST k-u/k-e model 

• Jones-Launder k - e model 

• Launder- Sharma k - e model (baseline model) 

• Chien k - e model 

• Baldwin-Lomax algebraic mixing-length model 

Four important model modifications are optional in the baseline 
models: 

• A length scale modification in order to decrease the heat transfer 
rate in flow reattachment zones (all models). 

• A rapid compression modification to increase the size of the flow 
separation zone in high speed flows (all models). 

• A rotation or vortex stretching length scale modification to de- 
crease the level of turbulent dissipation in flow recirculating 
zones (all models). 

• A turbulence transition modification to simulate flow transition 
and complex surfaces with both laminar and turbulent compress- 
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ible flow zones (k -u> model). 

The validation of the different models and modifications have been 
accomplish with comparisons against theoretical correlations, other nu- 
merical codes, and different experimental data. Flows included in the 
validation are: 

• Mach 14 laminar boundary layer on a 24 ° compression ramp 
(code validation). 

• Mach 5 flat plate turbulent boundary layer (turbulence modeling 
validation). 

• Mach 3 separated turbulent boundary layer on a 24 ° compression 
ramp (turbulent flow separation validation). 

• Mach 8 shock wave / turbulent boundary layer interaction gen- 
erated with a 15 ° and a 10 ° fin mounted on a cooled flat plate 
with Tuaii/Too = 0.3 and Re ^ = 5 • 10 6 (complex 3-D SWBLI with 
crossflow separation). 

• Mach 8 and 4 intersecting shock waves / turbulent boundary 
layer interaction generated with two 15 ° fins mounted on a flat 
plate (complex 3-D SWBLI with completely crossflow separa- 
tion). 

2. AIAA 94-1905 Paper (see appendix). 

“Three-Dimensional Navier-Stokes Simulations with Two-Equation 
Turbulence Models of Intersecting Shock- Waves/ Turbulent Bound- 
ary layer at Mach 8.3,” has been written and presented in the 12th 
AIAA Applied Aerodynamics Conference held in June 20-23, 1994, 
at Colorado Springs, Colorado. This paper is co-authored with Dr. 
T.J. Coakley. 
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This publication shows the validation and comparison of the baseline 
k-u two-equation model with and without different combinations of 
model corrections, the SST two-equation model, and the Baldwin- 
Lomax algebraic mixing-length model. It demonstrate the ability to 
predict accurately surface pressure, heat transfer rates, flow yaw an- 
gles and Pitot pressures. The model corrections give improved heat 
transfer predictions in the reattachment zones. A very important 
result is the influence of the fin boundary layer on the flat plate dis- 
tributions. The flow structures are well capture and superior results. 

3. AIAA 94-2950 Paper (see appendix). 

“Three-Dimensional Navier-Stokes Method with Two-Equation Tur- 
bulence Models for Efficient Numerical Simulation of Hypersonic 
Flows,” has been written and presented in the 30th AIAA/ ASME/ 
SAE/ ASEE Joint Propulsion Conference and Exhibit held in June 
27-29, 1994, at Indianapolis, Indiana. 

This publication shows the 3-D method used to predict the turbulent 
flows, and the necessary information to implement the two-equation 
turbulence models into 3-D CFD codes. The method is very efficient 
and allows the simulation of different turbulence models in complex 
3-D flows in less computational time than other well-known methods. 
This research has shown that the presence of oblique shock waves 
introduce errors in the upwind methods that propagate downstream 
and have a significant effects in the performance of the model equa- 
tions. A limiter has been designed and tested to avoid this numerical 
error. The use of implicit and coupled treatment of the model equa- 
tions provide a fast convergence in complex flows, including separa- 
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tion and reattachment of boundary layers. This paper provides the 
basic documentation and delivers the information on their numerical 
procedures needed to implement the baseline turbulence models into 
3-D CFD codes. 

4. NASA TM 108827 (see appendix). 

The final report of the GWP-18 in the Modeling and Experimental 
Validation Branch was written together with Tom Coakley. This 
report has been published as NASA TM 108827, “Turbulence Com- 
pressibility Corrections,” and includes the database, recomended 
models, and model equations for hypersonic flows. 

5. AIAA 95-2215 Paper (to be published). 

“The Structure of Intersecting Shock- Waves/ Turbulent Boundary 
layer Flow,” will be presented in the 26th AIAA Fluid Dynamics 
Conference, to be held June 19-22, 1995, in San Diego, California. 
This paper is co-authored with Dr. T.J. Coakley. 

This paper study the interactions of the shock- waves and the turbu- 
lence structures and the effects of the turbulence models at Mach 8 
in high resolution numerical simulations with over 1.5 million grid- 
points. This structures are currently not well understood, and their 
modeling requires a higher-level knowledgement of their physics. 
This simulations provides the best description of this flow and a very 
valuable database. Researchers in this area have already expressed 
their interest in this data set. The implicit boundary approxima- 
tions of the models have been tested and improved. The length-scale 
model correction show improvements in the prediction of peak sur- 
face pressure and heat transfer rate. The rapid-compression model 
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correction show increased secondary separation under the main vor- 
tical structures. These simulations study the effects and capabilities 
of different two-equation turbulence models in their predictions of 
complex turbulence. 

6. AIAA 95 Paper (to be published). 

The study of the validation of the different turbulence models with 
the intersecting shock- waves/turbulent boundary layer interaction 
flow at Mach 4 will be presented in an AIAA meeting to be held 
later this year. The organization of this meeting begun in 1994, the 
session chairman is J. Marvin and the final information will be soon 
available. This experimental database of the Garrison and Settles 
experiment provides surface presure and skin friction distributions, 
and flowfield Pitot pressure. The numerical simulations have been 
carried out with a 91x61x61 3-D grid and the baseline two-equation 
turbulence model with and without corrections. The few other sim- 
ulations available have consistently shown poor agreement with the 
experimental data in the skin friction coefficient (and heat trans- 
fer rate). The models also show significant differences in the flow 
streamline patterns. This meeting provides the mean to effective 
validation of models and methods in a complex turbulent flow. 

7. As part of a collaborative effort, several tests of the shock- tunnel inlet 
simulation have been run. Velocity profiles, temperature profiles, 
skin friction and heat transfer coefficient have been compared under 
different Reynolds numbers, Mach numbers, ratio of specific heat- 
transfer coefficients, and Prandtl numbers. 
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8. NASA Peer Review 

A Peer Review of Turbulence Modeling work at NASA was held in 
September. The review team consisted of Managers of aerospace 
companies and university Professors. The presentations showed the 
research work in progress at Ames (including CTR), Lewis, and Lan- 
gley. A set of slides was prepared for this meeting with compar- 
isons between models and experiment in 3-D shock- wave/ turbulent 
boundary layer interactions (flat plate/ 10° single fin, flat plate/ 15° 
single fin, and flat plate/ 15° double fin). The results show the effects 
of turbulence compressibility corrections in two-equation turbulence 
models in the prediction of surface pressure, heat transfer, skin fric- 
tion, yaw angles, and Pitot pressure measurements. This work was 
presented by Joseph Marvin. The review team found the research 
being excellent. However, they found a lack of general management 
of the whole NASA program. Their main recommendation was that 
turbulence research done at NASA is a fundamental need of the 
American industries, otherwise, the future technology will be lost to 
European or Japanese companies. Their recommendations include 
the appointment of a general manager of the turbulence and tran- 
sition programs with authority and outstanding expertize; creation 
of a board with members of industry and universities to guide this 
program; creation of a technical board to supervise the research; and 
a general coordination of all these NASA Center research activities. 
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Summary 


The basic objective of this research was to identify, develop and recommend 
turbulence models which could be incorporated into CFD codes used in the design of the 
NASP vehicles. To accomplish this goal, a combined effort consisting of experimental and 
theoretical phases was undertaken. The experimental phase consisted of a literature survey 
to collect and assess a database of well documented experimental flows, with emphasis on 
high speed or hypersonic flows, which could be used to validate turbulence models. Since 
it was anticipated that this database would be incomplete and would need supplementing, 
additional experiments in the NASA Ames 3.5' hypersonic wind tunnel (HWT) were also 
undertaken. The theoretical phase consisted of identifying promising turbulence models 
through applications to simple flows, and then investigating more promising models in 
applications to complex flows. The complex flows were selected from the database devel- 
oped in the first phase of the study. For these flows it was anticipated that model perfor- 
mance would not be entirely satisfactory so that model improvements or corrections 
would be required. The primary goals of the investigation were essentially achieved. A 
large database of flows was collected and assessed, a number of additional hypersonic 
experiments were conducted in the Ames HWT, and two turbulence models (k-e and k-co 
models with corrections) were determined which gave superior performances for most of 
the flows studied and are now recommended for NASP applications. 

Introduction 

With the advent of the hypersonic airplane, hypersonic flows are receiving spe- 
cial attention from researchers in computational fluid dynamics. Complex flow phenom- 
ena such as shock-wave boundary layer interactions, separation and combustion are of 
particular interest because of their importance to the successful design of structural, pro- 
pulsive and thermal protection systems. Rapid advances in CFD in recent years have 
resulted in its increased use as a design tool for aeronautical systems and have lead to 
reductions in the time and costs of wind tunnel testing. A major obstacle to the use of CFD 
as a design tool is its dependence on turbulence modeling for accurate prediction of com- 
plex flows. Although recent advances have been made in turbulence modeling, many more 
will be required before CFD can be applied with confidence to a wide range of flow prob- 
lems. This is especially true at hypersonic speeds where high temperatures and pressures 
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create additional difficulties for turbulence modeling and where wind tunnel experiments 
which can be used to validate models are sparse. 

The overall goal of this research is to identify, develop and recommend turbu- 
lence models which can be incorporated into the CFD codes used in the design of the 
NASP vehicle. To accomplish this objective, a combined and joint effort consisting of 
experimental and theoretical phases was undertaken. The objective of the experimental 
phase was to conduct a literature survey to identify and assess a database of well docu- 
mented experimental flows which could be used to validate turbulence models. Since it 
was anticipated that this database would be incomplete, it was also decided to perform 
additional experiments in the NASA Ames 3.5’ hypersonic wind tunnel (HWT) to supple- 
ment the database. 

The objective of the theoretical phase was to identify promising turbulence mod- 
els through applications to simple flows, such as flat plate flows, and to then investigate 
the more promising models in applications to more complex flows. The complex flows 
were to be selected from the database developed in the first phase of the study. For these 
flows it was anticipated that model performance would not be entirely satisfactory so that 
model improvements or corrections would be required. The flows of interest were 
restricted to ideal gas flows because of the sparsity of high quality experimental validation 
data and viable turbulence models for real gas flows. 

The schedule of tasks and milestones for the completion of the research on Gov- 
ernment Work Package 18 is shown in fig. 1. With this report the work is essentially com- 
plete with the exception of the compressible shear layer experiment in the Ames 3.5’ 
HWT. Completion of this experiment was halted due to the lack of funds caused by fund- 
ing reductions in the NASP project. 

The report is organized into 7 sections. Following the introduction, two sections 
on the experimental phase of the study will be presented including one on the database 
collection and assessment and another on the experiments conducted in the Ames 3.5’ 
HWT. Next, sections describing the theoretical phase of the study will be presented. These 
include sections on recommended baseline turbulence models, compressibility corrections 
recommended for improved predictions of complex flows and representative results of 
numerical predictions using the baseline and corrected models. Finally, the report con- 
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eludes with a summary of basic results and recommendations including topics for future 
study. 


Database Collection and Assessment 

As stated in the introduction, the purpose of the database collection and assess- 
ment activity was to provide a base of reliable and well documented wind tunnel experi- 
ments by means of which turbulence models could be validated. These experiments, for 
the most part, involve relatively simple geometric shapes which may be viewed as sepa- 
rate elements of an overall vehicle. The results of this activity are reported by Settles and 
Dodson (1991, 1993a and 1993b) and will be summarized here. 

Settles and Dodson (1991 and 1993b) provide a survey and assessment of exper- 
iments involving two and three dimensional shock- wave boundary-layer interaction flows. 
Eight hundred experiments were initially identified for further review. Of these, 112 dis- 
tinct experiments were found involving flows at Mach 3 and above. The acceptance crite- 
ria applied to these included: 1) Measurements of surface pressures, skin friction and/or 
heat transfer, and velocity, temperature, or pitot pressure profiles at selected locations, 2) 
well defined experimental boundary conditions, 3) well defined error bounds, 4) adequate 
spatial resolution of measurements and 5) full documentation of tabulated data. For hyper- 
sonic conditions, i.e. M > 5, only 7 studies passed the acceptance criteria and only three 
were three dimensional. An additional 11 experiments passed the criteria at supersonic 
speeds. 


The survey and assessment of attached boundary layer and free shear flow exper- 
iments is given by Settles and Dodson (1993a). The acceptance criteria applied to these 
cases was identical to that applied to the shock-wave boundary-layer interaction experi- 
ments. For the boundary layer experiments, 153 candidate cases were identified for further 
review. Of these 39 were subjected to the acceptance criteria. No hypersonic and only 9 
supersonic cases passed. For the free shear layer experiments, 1137 candidate cases were 
identified for further review. Of these, 45 were subjected to the acceptance criteria and 
only 3 passed. 
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Experiments in the NASA Ames 3.5’ Hypersonic Wind 

Ttinnel 

The experiments were conducted in the Ames 3.5’ hypersonic wind tunnel and 
were done with relatively simple generic shapes such as cones, cylinders, plates and 
wedges. These shapes were chosen to typify locations on a high speed vehicle where tur- 
bulence modeling was expected to be a critical issue. The experiments were run at nomi- 
nal Mach numbers ranging from 7 to 8.3 and unit Reynolds numbers (per meter) from 4.9 
to 5.8 million. Boundary-layers approaching the interaction region were relatively large, 
with thicknesses on the order of 2.5 to 3.7 cm, which allowed detailed flow field surveys to 
be easily made. In most of the experiments, both surface measurements and flow field 
(profile) measurements were made including initial boundary conditions required to start 
numerical computations. With two minor exceptions, the experiments passed the accep- 
tance criteria described in the previous section. In addition, an analysis of measurement 
errors was made and documented. The results of these experiments are described in detail 
by Kussoy et al (1989; 1991a,b; 1992; 1993a, b,c) and Horstman and Kussoy (1989) and 
the data have been made available on floppy disks. 

M = 7 Hypersonic Cylinder-Flare and Fin flows 

The test bed employed in this experiment consisted of an ogive-cylinder at zero 
angle of attack with a series of removable symmetric flares or sharp fins (see figs. (2a and 
2b)). Both flare and fin angles were varied, producing shock waves of various strengths, 
and resulting in both attached and separated flow fields. Detailed measurements verified a 
fully developed turbulent boundary- layer on the cylinder ahead of the interaction region. 
The resulting flows were axisymmetric with and without separation for the flare case, and 
three dimensional with separation for the fin case. Surface pressures and heat transfer rates 
were measured on both configurations, and flow field surveys were done on the flare con- 
figuration. The results are reported in Kussoy and Horstman (1989) and Horstman and 
Kussoy (1989). 


M = 8.2 2-D Wedge and 3-D Vertical Fin Flows 

A flat plate arrangement was used for this experimental series and is shown in 
fig.(3). it was of a hollow modular construction, enabling both test bodies and instrumen- 
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tation to be easily manipulated and changed. The full length of the wind tunnel test section 
was utilized in the design and because of this a well developed equilibrium turbulent 
boundary-layer was present at the shock interaction zone. Two configurations were tested; 
the first consisted of a sharp wedge supported over the width of the test section, the second 
was a sharp vertical fin attached to the plate surface. These are both illustrated in fig.(3). 
Both the wedge and fin angles were varied, producing shock-wave boundary-layer interac- 
tions of varying strength with a maximum wall-pressure ratio of pjp^ = 21.5 and 6.4, 
respectively. This resulted in both attached and separated flow fields for the wedge flows, 
and swept three dimensional vortical flow fields for the fin flows. Detailed surveys verified 
a fully developed hypersonic turbulent boundary-layer on the flat plate ahead of the inter- 
action zone. 

For the wedge configuration, only surface conditions were measured. For the fin 
configuration, however, mean flow profile surveys were also taken - both in the undis- 
turbed and interaction regions - and from them pitot pressure contours and boundary-layer 
thickness parameters were obtained. We believe this is the first fully three dimensional 
shock-wave boundary-layer interaction flow to be so documented at hypersonic speeds. 
Results for both configurations are reported in Kussoy et al (1991a,b;1992). Experimental 
and computational results for the 10° and 15° vertical fin flows are discussed in the sec- 
tion on Model Validation-Representative Results. 

M = 8.3 Crossing Shock Flow 

For the third series of experiments, a configuration was chosen to reflect several 
key elements of a generic hypersonic inlet. These included a thick turbulent boundary- 
layer approaching two vertical fins of varying wedge angle, a crossing shock pattern pro- 
duced by the fins, boundary-layer vortices, large pressure gradients, and separation zones. 
The test body for this series of experiments is shown in fig.(4). 

Streamwise and transverse surface pressure and heat transfer distributions were 
measured as well as flow field surveys of pitot pressure and flow angle. One important 
result of these measurements should be mentioned here. This was the persistence of an 
extensive low pressure region far downstream of the fin leading edges. This low pressure 
region implied that the generic inlet tested here would not be a very efficient pressure dif- 
fusing device. The experimental results for this configuration are given in Kussoy et al 
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(1993a,b,c). experimental and computational results for the 15° fin-angle case are dis- 
cussed in the section on Model Validation-Representative Results. 

Recommended Baseline Turbulence Models 

The baseline turbulence models recommended for use in NASP applications will 
be described in this section. A variety of turbulence models have been investigated in 
varying degrees throughout the course of the study. They include O-eq, 1-eq, and 2-eq 
eddy viscosity models, and Reynolds stress transport models. For NASP applications, the 
primary emphasis has been placed on 2-eq models and these models will be the only ones 
described in detail here. Descriptions of other models used in the course of the study, and 
results obtained with them, are given in Coakley and Huang (1992) and Coakley and Mar- 
vin (1993), Hors tm an (1991,1992), and Huang and Coakley (1993a, b). 

The models investigated and recommended here are those that utilize no slip 
boundary conditions at solid walls and involve the use of wall damping functions (in most 
cases). This is in contrast to conventional practice with two equation models in which wall 
functions and essentially slip type boundary conditions are used. This choice was made 
because at the very high speeds of hypersonic flight, the effective Reynolds numbers of 
the flows can be quite low, in some cases involving transition and relaminaiization. In 
these cases, the thickness of the laminar sublayer becomes an appreciable fraction of the 
overall boundary layer thickness and the wall function approach becomes inapplicable or 
ineffective. In addition, the wall function approach gives questionable results for separated 
flows. For these reasons it was decided to use the wall damping function approach. 

Two equation models have been emphasized since these are viewed as the sim- 
plest and most practical models available which have sufficient generality to be applied to 
the complex flows of interest in NASP applications. Although numerous two equation 
models have been investigated, descriptions of only two of these will be given here since 
these are the baseline models of our final recommendation. In Coakley and Huang (1992) 
a detailed investigation of the performance of a variety two equation models was pre- 
sented for flat plate boundary layers over a wide range of Mach and Reynolds numbers. 
The results of that study showed that most of the models gave reasonably good predictions 
of skin friction, heat transfer and velocity profiles with little clear preference of one model 
over the other. For this reason only two models were selected for further study, and it is 
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believed that these models are representative of most two equation models currently in 
use. The baseline models are the k-e model of Jones and Launder (1972), as modified by 
Launder and Sharma (1973), and the k-a> model of Wilcox (1984). These models are used 
with the (mass weighted) Reynolds averaged compressible Navier-Stokes equations 
which, in cartesian tensor form, are given below. 

Reynolds Averaged Compressible Navier-Stokes Equations 

|(p) + A(p« > )=° 

+ 57 ( P",“j + °,y> =0 

S <p£) + S / (p£ V“.<V'?, ) =0 


where; p is the density; w, are the cartesian velocity components; E = e + + k is the 

total specific energy; e = c v T is the specific internal energy; k = 05 pufuf /p is the turbu- 
lent kinetic energy; T is the temperature; p = (y - l)pe is the equation of state; and p is 
the pressure; y = c pi c v is the ratio of specific heats, and c p and c v are the specific heats at 
constant pressure and volume respectively. The variables Gjj and qj are the total stress 
tensor and heat flux vector, respectively, which include both molecular and (Reynolds 
averaged) turbulent contributions. Using the Boussinesq approximation, these variables 
are represented in terms of an eddy viscosity by 


Oij = h i j(p+^pk)-(\L + [L T ) 

r U M BT 

^CT Oj) 


fdu t Buj 
ydXj + Bx 


<7 . = - y 


H‘5) 

;-R)l 


where |! and |i r are the molecular and turbulent (eddy) viscosities a and oj are molecular 
and turbulent Prandtl numbers with o = c p JJ/k (assuming air) o T = 0.9 and depending 

on the model used. The turbulent eddy viscosity is expressed in terms of the turbulent 
kinetic energy, k, and either the dissipation rate, e, or the specific dissipation rate go 
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depending on the model. This expression is 


r k 1 k e 

\kj — Cyf jkl = Cyfyfi— — Cyfy jP— , (o — 

where V* is the turbulent velocity scale, / = V/rfye = V*/© is the length scale, is a mod- 
eling constant, and is a damping function depending on the specific model used. For 
applications to complex flows, the governing equations were expressed in terms of curvi- 
linear coordinates and solved using the finite volume method, Viegas and Rubesin (1991), 
and Huang and Coakley (1992b) or the finite difference method, Bardina (1994). 

The recommended baseline k-e and k-a> turbulence models are expressed by the 
formulas given below. 

h **> 

a ( PJ ) + £-<p su r ^) . (P,( S -) 2 -a?-D,)pw 
_ /dUf duftdui- 2(<hk) 2 n 

^dxj + 3xj5xJ 3^dxJ ’ dx k 

In these equations and in the following Tables, variable “s” and subscript “s” are replaced 
with e for the k-e model or © for the k-ca model, respectively. 

The model parameters are defined in Tables 1, 2 and 3 


Table 1: Model Parameters 


= 9/100 

\l k = \L + \Lrlo k 

\L s = \l + \lrlo s 

R t = k?/ve = k/v(o 


II 

£ 

a s = 2/3 C sl 

Cl 

II 

> 

v T = \h' l P 
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Table 2: Launder-Sharma k-e Model, s ■ e = coyt 


1 — — 

CT*=1 

1-3 

C El = 1.45 

Cg2 = 1.92 

fp = exp( - 3.4/(1.0 + R t 2 /50)) 

h - 1 - 0.3 exp( - R t 2 ) 

D k = 1 + (2v/e) (dk JI2 /dx k ) 2 

D e = 02 / 2 - (2v Vj/e 2 ) (<P(u k u k ) ll2 /dx k dx k ) 2 


Table 3: Wilcox k - co Model, 5 s © = ejk 


$ 

It 

^3 

^=1/2 

c al = 5/9 



ESI 

3 

II 

3 

Q 


For the Launder-Sharma model, the thin layer approximation was normally used 
to compute the derivatives in the D k and D e terms. 

An important consideration in using the Wilcox k-ca model, which is not neces- 
sary with the k-e model, is the value of to in the free stream (just outside the boundary- 
layer edge) where it cannot be too small. In all of the applications of this study, it was pos- 
sible to choose values of co at the inflow boundary which insured that the values of co in 
the free stream would not be too small. (It must be chosen such that C 0 oo> 10 UJL where 
Uoo is the free stream velocity and L is the length of run of the boundary layer and cOoo is 
the free stream value of co at the start of the boundary layer). It may be that in future appli- 
cations it will not be feasable to control the free stream 00 oo in this manner, and other mea- 
sures will be necessary. One alternative would be to use the Jt-co model of F. Menter 
(1992), which uses a blending of the k-ca and k-e models to circumvent the problem. This 
model will be discussed more fully in the next section. 

Recommended Model Corrections for Complex Flows 

The baseline models described above are generally not adequate to accurately 
predict complex flow problems and must be corrected to deal with these cases. The types 
of modeling corrections that have been found useful in practice and which are recom- 
mended for NASP applications are summarized below. Other model corrections which 
have been tried in the course of the study but which are not recommended for applications 
are described in Coakley and Huang (1992). 


9 















Length Scale Correction 


The first correction is addressed to difficulties encountered in predicting heat 
transfer in the reattachment or shock impingement zone of shock-wave boundary-layer 
interaction flows. In these zones, all 2-eq models dramatically overpredict heat transfer 
and must be corrected. The correction involves the use of an algebraic length scale which 
limits the slope of the length scale predicted by the two equation model, which otherwise 
would become very large in these regions. The formulas defining the length scale correc- 
tion for the k-e and £-© models are given below, 

/ = min{2.5y,J^/e) = min{2.5y,Jk/(i)} 

in these formulas, / is the turbulent length scale which is taken to be the smaller of an alge- 
braic expression (kC ^' 3/4 y = 2.5 y, based on a von Karin an constant of k = 0.41) and the 
conventional length scale given by the two equation model. Having computed this param- 
eter, the value of e or © is recomputed and reset to be consistent with this value, e.g. e = 

V* 3 // or o>= V*//. 


Rapid Compression Correction 

The second correction to be described is called the rapid compression correction 
and is used to improve predictions of separation in shock-wave turbulent boundary-layer 
interactions. The correction involves changing the coefficient of the dilitation or velocity 
divergence in the e and © equations, i.e. ot^ or a^, The net effect of this correction is to 
increase the production of epsilon or omega in regions of rapid compression, or shock in 
waves, which reduces the eddy viscosity and enhances separation. The corrected values of 
the dilitation coefficient for each model are 

a £ = 2 or a ffl = 4/3 

as opposed to the Launder-Sharma coefficient a £ = (2/3) C^j - 0.97 and the Wilcox coeffi- 
cient a w = (2/3) C m j = 0.37 shown in Tables 1, 2, and 3. The development of this correc- 
tion is discussed more fully in Coakley and Huang (1992). 
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Compressible Dissipation Correction 


The third and final correction recommended for NASP applications is applicable 
to free shear flows (e.g. mixing layers, wakes and jets). In these types of flows, it is well 
known that shear layer spreading rates decrease as Mach numbers increase compared with 
spreading rates at zero Mach number. The correction recommended to improve predic- 
tions of these flows was developed originally by Zeman (1990), and is closely related to 
similar corrections developed by Sarkar (1991), and Wilcox (1992). The correction is 
listed below. 


D k^ D k + $k 
D a> ~ $k 

= M 1 - ex P (-(max (0, a x M T - a 2 ) /a 3 ) 2 ) ) 

D e is unchanged, M T = V*/c, c is the local sound speed, a 0 = 3/4, aj = V(yf 1), a 2 = 1/10, 
and aj = 6/10. 

For applications involving boundary layers, or solid walls, this correction has 
been found to underpredict skin friction especially at high free stream Mach numbers. It is 
recommended, therefore, that this correction only be used in free shear flow applications. 
The correction of Wilcox (1992), which is similar to that of Zeman, has been found to 
work well in boundary layers as well as free shear flows and Wilcox recommends its use 
without reservation. Since we have not investigated this model under the wide range of 
conditions investigated using other model corrections, we chose not to recommend the 
model at this time. 


Other Corrections and Models of Interest 

Although the above models and corrections constitute our final recommenda- 
tions for NASP applications, it is important to note that these models only constitute an 
improvement over previous models and may very well give poor predictions for flow situ- 
ations which have not yet been investigated and validated. In this regard we draw attention 
to the k - co model of F. Menter (1993), which has proven quite successful in incompress- 
ible and transonic flow applications. This model was designed to overcome certain defi- 


ll 



ciencies of the Wilcox k-a> model and has certain features which might prove useful in 
hypersonic applications. One feature of the model is the so called shear stress transport or 
rapid strain feature which enables the model to give improved predictions of adverse pres- 
sure gradient boundary-layers and separation. This feature is similar to the rapid compres- 
sion correction discussed above and enhances separation by reducing the eddy viscosity in 
non-equilibrium regions where the flow is changing rapidly. This feature was also tried 
with the baseline models and gave results very similar to those obtained with the rapid 
compression correction. It is believed that some combination of the two corrections may 
ultimately prove more accurate and reliable in future applications. It must be stated, how- 
ever, that we did apply Menter’s model to most of the flows described in this study and 
obtained no improvement over the recommended models. In some cases it did not perform 
as well. Since the model is considerably more complicated than the other models studied 
we decided not to include it in the list of recommended models. 

Model Validation - Representative Results 

Representative results of calculations and comparisons of model predictions with 
experimental measurements will be given in this section. The flows discussed include free 
shear (mixing layer) flows, 3 two dimensional shock-wave boundary- layer interaction 
flows, and 2 three dimensional shock-wave boundary-layer interactions. The turbulence 
models used include the recommended baseline and corrected k-E and k- co models. 

Compressible Mixing layer 

The first series of flows to be discussed consists of high speed mixing layers 
which are of considerable importance in the design of propulsive exhaust nozzles for the 
NASP vehicle. The comparisons of computations with experimental measurements is 
shown in fig.(5). The figure shows predictions of spreading rate divided by spreading rate 
at zero Mach n um ber compared with experimental measurements over a range of convec- 
tive Mach numbers. The data include the Bogdanoff (1993) compilation of the Langley 
fin ta , and the measurements of Samimy and Elliot (1990). The calculations were done by 
Viegas and Rubesin (1992), who used the baseline k-e model with the compressible dissi- 
pation corrections of Zeman and Sarkar. Although not shown, results obtained with the k- 
(0 model give similar results. Examination of the comparisons indicates that the baseline 
model significantly overpredicts the spreading rate while the model corrections improve 


12 



the predictions. Of the two model corrections, the Zeman correction gives the best overall 
result. 

Two-Dimensional Shock-Wave Boundary-Layer Interaction Flows 

The second series of model validation studies to be discussed consists of two 
dimensional (planar and axisymmetric) shock-wave boundary-layer interactions. The 
flows are the Mach 7 ogive-cylinder-flare flow of Kussoy and Horstman (1989) (35° flare 
angle), the Mach 9 planar compression ramp flow of Coleman and Stollery (1972) (34° 
ramp angle), and the Mach 7 axisymmetric impinging shock flow of Kussoy and Horst- 
man (1975) (15° generator angle). In all cases, the walls were highly cooled with wall-to- 
adiabatic wall temperature ratios on the order of 0.3 to 0.4. The test configurations for 
these cases are shown in figs.(6a,b,c). 

All calculations were done using the code developed by Huang and Coakley 
(1992). The inlet flow conditions just ahead of the shock interaction zones were obtained 
by calculating the flow over a flat plate and matching measured and computed displace- 
ment thicknesses. The value of y + at the first grid point off the wall was maintained to be 
less than 0.5 and the grid was expanded exponentially from the wall to the free stream. 
This gave between 60 to 80 grid cells in the boundary-layer and 140 cells overall in the 
cross stream direction. Computations with fewer cells inside the boundary-layer (i.e. 40 
cells) were made and no significant differences were observed. In the stream wise direction 
a grid of 140 cells was used except in the impinging shock case where 200 cells were used. 
The comparisons of computations with measurements are shown in figs. (7-9). They 
include measured and computed surface pressures and heat transfer distributions (and skin 
friction for the impinging shock case) and were done with the baseline and corrected k-e 
and k-(o models. The pressure and heat transfer measurements are normalized by the mea- 
sured values in the region ahead of the interaction. Zeman’s compressible dissipation cor- 
rection designed for shear layers was not used as explained in the section on model 
corrections. It is clear from these results that the baseline models significantly under pre- 
dict the extent of separation and over predict the heat transfer in the interaction region. 
The corrected models both give results in much better agreement with experiment. 
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Three-Dimensional Shock-Wave Boundary-Layer Interaction flows 


The Ames experiments on 3-D shock-wave boundary-layer interactions 
described in this report were used here to test the compressible turbulence models and the 
model corrections. The experiments and the data are also described in Settles and Dodson 
[1993a and 1993b] and Kussoy and Horstman [1991,1992,1993]. The calculations were 
done with the code of Bardina (1994) 

3-D Vertical Fin Shock Interaction Flows 

This experiment investigates the interaction of a hypersonic shock wave with a 
thick turbulent boundary layer [Kussoy and Horstman, 1991 and 1993b]. A 10° and a 15° 
vertical fins mounted on top of a flat plate were used to generate oblique shock waves. The 
free-stream Mach number was Mo = 8.2, the temperature was To = 81° K, and the Rey- 
nolds number was Re« = 5’ 10 6 per meter. The wall temperature was fixed at 300° K. The 
interaction of the shock wave with the turbulent boundary layer generates a crossflow vor- 
tex separation with a “quasi conical’’ shape [Settles and Lu, 1985; Knight, Horstman, and 
Monson, 1992], Peak wall pressure, skin friction, and heat transfer rates were observed in 
the re-attachment zone behind the crossflow vortex. 

The numerical computations were made with 61x41x61 an 31x21x31 meshes 
[Bardina, Coakley, and Marvin, 1992]. Only small differences between the solutions were 
observed, and the fine mesh solutions are considered accurate for engineering purposes. 
The inflow conditions were obtained from the Navier-Stokes code solution matching the 
experimental displacement thickness. 

A few comparisons of experiment and simulation with the finer mesh are 
described below. Figures 10a and 10b show the surface pressure and skin friction distribu- 
tions, respectively, for the 10° fin flow on the flat plate surface at the crossed section 
located at x=0.1819 m downstream of the fin leading edge. Figure 10c shows the compar- 
ison of the wall heat transfer distribution on the flat surface at x=0.1645 m downstream of 
the fin leading edge. Comparable results for the 15° fin case are shown in fig.(lla,b,c). 
The symbols in the figures show the experimental data points, the solid lines show the 
solution with the baseline £-co model, and the dash lines show the solution with the k- to 
model with both model corrections (length-scale and rapid compression corrections). Both 
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simulations show good agreement with the experimental data. The peak values observed 
in the re-attachment zone are also well predicted. The simulations fail to predict the small 
plateau observed in the wall heat transfer rate distributions upstream of the shock wave. 
Comparison of baseline and corrected model predictions show only small differences. 
This is probably due to the fact that, compared with the 2-D results, the pressure rise 
through the shock wave is relatively weak and separation is relieved by three dimensional 
effects. 


3-D Crossing Shock Interaction Flow 

This experiment studies the interactions of two intersecting hypersonic shock 
waves with a thick turbulent boundary layer [Kussoy and Horstman, 1992]. Two 15° fins 
mounted on top of a flat plate were used to generate intersecting oblique shock waves (see 
fig. 4). The free-stream Mach number was Moo = 8.3, the temperature was Too = 80° K, and 
the Reynolds number was Reoo = 5.3T0 6 per meter. The wall temperature was fixed at 
300° K. The intersection of both crossflow vortices generates different complex flow 
structures with high static pressures and surface heat transfer rates. The intersection of the 
two “quasi conical” vortical structures uplifts the flow and induces a wave structure in the 
symmetry plane [Gaitonde and Shang, 1993]. 

The numerical simulations were made with a 231x81x81 mesh [Bardina and 
Coakley, 1994; Bardina, 1994]. Simulations studies with 101x61x41 and 31x21x31 grid 
points were also done to analyze grid effects. Small differences between the solutions 
were observed in the surface pressure and heat transfer distributions, but signifficant dif- 
ferences were observed in flow structure. The fine mesh solutions provided the best reso- 
lution of the turbulence structures and are considered accurate for engineering purposes. 
The inflow conditions were obtained from the Navier-Stokes code solution matching the 
experimental displacement thickness. 

Figures 12a and 12b show the pressure and heat transfer distributions, respec- 
tively, on the plate surface along the symmetry plane located between the two fins. The 
symbols show the experimental data, the solid lines show the solution with the baseline k- 
co model, and the dash lines show the solution with the Jt-co model with both model correc- 
tions (length-scale and rapid compression corrections). Predictions of both surface pres- 
sure and heat transfer show good agreement, within the experimental uncertainty, except 
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near the outflow zone. The small plateau at the beginning of the interaction is also not pre- 
dicted by either model. The peak pressure and heat transfer rate are very well predicted. 
As in the single fin cases, both the baseline and corrected model predictions show only 
small differences. 

Figures 13a and 13b show the pressure distributions on the plate surface at two 
cross-sections, one located at x/boo = 5.60 and the other located at x/5oo = 6.92 downstream 
of the fin leading edge. The first distribution is upstream of the peak surface pressure, and 
the second distribution is downstream of the peak surface pressure generated by the flow 
re-attachment. Both simulations show good agreement within the experimental uncer- 
tainty of the data, and both models give similar predictions. 

Figures 13c and 13d show the comparison of heat transfer rate profiles on the 
plate surface in the crossed sections located at x/Soo = 5.08 and x/5oo = 6.40, respectively, 
upstream and downstream of the peak heat transfer rate observed in the symmetry plane. 
Both simulations agree in general with the experimental data. The model corrections show 
improved heat transfer rate predictions in the downstream zone. 

Figure 14 shows experimental pitot pressure contours compared with computa- 
tional contours obtained with the corrected k - © model. Three locations are shown, x/5oo = 
5.60, 6.9, and 8.3, respectively. The resolution of the numerical data is 81x81 grid points 
while the resolution of the experimental data is 5x24, 4x24, and 4x24 respectively. The 
agreement with the experimental data is very good. 

Summary and Conclusions 

In this section we summarize the work performed under GWP 18, give our prin- 
cipal results and recommendations, and discuss plans for future work. We feel that, over- 
all, the results produced in the course of the work were of a very high caliber and will be 
of considerable use to modelers and designers of hypersonic flight vehicles. The research 
was divided into two distinct phases; one experimental and the other theoretical. The 
experimental phase consisted of the collection and assessment of a database of high speed 
wind tunnel experiments gathered from sources around the world and the conduct of addi- 
tional experiments in the NASA Ames 3.5’ hypersonic wind tunnel which would add to 
the database. The primary aim of this effort was to produce a database of reliable and rele- 
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vant flow measurements which could be used to validate turbulence models for NASP 
CFD design codes. Out of a total of over 2000 experiments that were initially studied, the 
list of acceptable experiments was reduced a total of 30 using a set of strict acceptance cri- 
teria. The kinds of flows analyzed and recommended as legitimate candidates for inclusion 
in the database included two- and three-dimensional shock-wave boundary-layer interac- 
tion flows, attached boundary layers flows and free shear flows. 

The theoretical phase of the work consisted of identifying, testing, and recom- 
mending turbulence models which would be of practical use in the CFD design codes. A 
large number of baseline models were initially tested on flat plate flows and this number 
was then reduced to two for further testing on more complex flows. The baseline models 
selected were the Launder-Sharma version of the k-e model, and the Wilcox k-i o model. 
These models were tested on free shear flows and two- and three-dimensional shock-wave 
boundary-layer interaction flows. Some of the results of these calculations have been dis- 
cussed in this report. It was found that the baseline models did not perform satisfactorily 
with regard to separation and heat transfer predictions especially for the 2-D shock-wave 
boundary-layer interaction flows, and did not accurately predict the spreading rate of free 
shear layers. To improve model performance for these complex flows, a series of com- 
pressibility corrections was investigated. The more promising of these corrections were 
then selected to be the final model recommendations for incorporation in to the NASP 
CFD codes. Representative results using the corrected models were discussed in this 
report. It was shown that for the flows investigated, the model corrections give substan- 
tially improved predictions. 

Although the work performed in the course of this research has led to the identi- 
fication of useful flows for model validation and the development of improved turbulence 
models for hypersonic flight, much work remains to be done. Not all of the flows included 
in the database of recommended flows have been investigated computationally, and these 
need to be investigated. More complicated flows not included in the database also need to 
be investigated. Flows that fall into this latter category include flows with chemical reac- 
tions such as those occurring in SCRAM jet combustion and propulsion, and flows that 
include transitional phenomena just to mention two. Interest in these and other areas of 
turbulence modeling research is very high and the work goes on. 
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Fig. 5. Effect of compressibility on normalized vorticity thickness growth rate for mix- 
ing layer (Viegas and Rubesin (1991, 1992)). 
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Fig. 6a. Axisymmetric compression generated with ogive-cylinder-flare configuration of 
Kussoy and Horstman (1989). 



Fig. 6b. 2-D planar compression ramp configuration of Coleman and Stollery 

(1972). 


Fig. 6c. Axisymmetric impinging shock generated with a cone-ogive cylinder and an 
wnnnlar fin configuration of Kussoy and Horstman (1975). 
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Fig. 7. Surface pressure (a) and heat transfer (b) in axisymmetric compression 
flow of Kussoy and Horstman (1989), geometry shown in Fig. 6a. 
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Fig. 8. Surface pressure (a) and heat transfer (b) in 2-D planar compression ramp 
flow of Coleman and Stollery (1972), geometry shown in Fig 6b. 
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Fig. 9. Surface pressure (a), heat transfer (b) and skin friction (c) in axisymmet- 
ric impinging shock flow of Kussoy and Horstman (1975), geometry shown in Fig 6c. 
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Fig. 11. Surface flow predictions in the 3-D vertical fin interaction flow (a=15° 
case) of Kussoy and Hors tm an (1991a), geometry shown in Fig. 3. 


33 







Fig. 12. Surface flow predictions in the symmetry plane of the 3-D crossing shock 
flow (a=15° case) of Kussoy and Horstman (1992), geometry shown in Fig. 4. (Symbols 
represent experimental data, — Jt-oo model, and — k - ca model with corrections). 
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Fig. 13. Surface flow quantities in the 3-D crossing shock interaction flow (a=15 
case) of Kussoy and Horstman (1992), geometry shown in Fig. 4. (Symbols represent 
experimental data, — k - to model, and — k - to model with corrections). 
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Fig. 14. Pitot pressure contours on cross-sections of 3-D crossing shock interaction 
flow (a=15° case) of Kussoy and Horstman (1992), geometry shown in Fig. 4. (The fin 
leading edges are located at x = 0, and 5oo is the thickness of the undisturbed turbulent 
boundary layer at x = 0). 
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Abstract 

An investigation of the numerical simulation with 
two-equation turbulence models of a three-dimensional 
hypersonic intersecting (SWTBL) shock-wave/turbulent 
boundary layer interaction flow is presented. The flows are 
solved with an efficient implicit upwind flux-difference 
split Reynolds-averaged Navier-Stokes code. Numerical 
results are compared with experimental data for a flow at 
Mach 8.28 and Reynolds number 5.3- 10 6 with crossing 
shock-waves and expansion fans generated by two lateral 
15 fins located on top of a cold-wall plate. This experiment 
belongs to the hypersonic database for modeling valida- 
tion. Simulations show the development of two primary 
counter-rotating cross-flow vortices and secondary turbu- 
lent structures under the main vortices and in each comer 
singularity inside the turbulent boundary layer. A signifi- 
cant loss of total pressure is produced by the complex 
interaction between the main vortices and the uplifted jet 
stream of the boundary layer. The overall agreement 
between computational and experimental data is generally 
good. The turbulence modeling corrections show improve- 
ments in the predictions of surface heat transfer distribu- 
tion and an increase in the strength of the cross-flow 
vortices. Accurate predictions of the outflow flowfield is 
found to require accurate modeling of the laminar/turbu- 
lent boundary layers on the fin walls. 

Introduction 

The present investigation is a continuous research 
effort to develop, verify and appily two-equation turbu- 
lence models for three-dimensional compressible turbulent 
flows 1 ’ 8 . Two-equation turbulence models are simple, gen- 
eral, robust, and efficient for CFD applications on high 
speed flows. Their simplicity is mainly due to the eddy 
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viscosity hypothesis; their generality is built on the use of 
transport equations to define turbulent velocity- and 
length-scale instead of algebraic definitions; their robust- 
ness is based on recent advances on numerical methods; 
and their efficiency is based on continuous improvements 
of present numerical processors. There is clear evidence 
showing that most widely used two-equation models tend 
to under-predict flow separation and over-predict heat 
transfer near flow re-attachment regions. In hypersonic 
flow calculations, these model deficiencies are even more 
pronounced, particularly regarding their inability to pre- 
dict the extent of the flow separation. 

In a recent investigation, Bardina and Coakley 1 tested 
two model corrections that were designed to remedy the 
above mentioned difficulties for 3-D hypersonic flows. 
Previously, Coakley and Huang 2 tested these model cor- 
rections against experimental data in 2-D compressible 
flows. The first one limits the turbulence length scale to be 
no greater than the von K£rm£n length scale. This correc- 
tion is equivalent to the use of a one-equation model in 
regions where the length scale of the two-equation model 
is larger than the von Karm&n length scale. The main 
effect of this correction was observed to reduce the heat 
transfer rate near flow-reattachment in agreement with 
experimental observations. The second model correction, 
designed to increase the extend of separation, causes the 
length scale to decrease (or increase) when the flow under- 
goes rapid compression (or expansion). 

The development and validation of turbulence models 
for hypersonic shock-waveAurbulent boundary layer inter- 
actions is based on the fundamental understanding of the 
turbulence physics of the flows and the availability of an 
acceptable experimental database. In a parallel research 
effort to this investigation. Settles and Dodson 9 ' 11 have 
completed an extensive review of the available experi- 
mental data on compressible turbulent high-speed flows 
suitable for turbulence model validation. This research 
was developed for the (NASP) National AeroSpace Plane 
program in the Modeling and Experimental Validation 
Branch at NASA Ames Research Center. The review of 2- 
D and 3-D flows includes complex hypersonic flows with 
pressure profiles, skin friction, wall heat transfer, yaw 


1 



angles, pitot pressure, and turbulence statistics data. In 
this paper, we present comparisons of numerical simula- 
tion results and the selected experimental data of Kussoy 
and Horstman 12 on intersecting shock- waves/ turbulent 
boundary layer (ISWTBL) interactions. This simple 
geometry shows complex turbulence structures of great 
interest in the design of uniform high-pressure flows at the 
entrance of inlets. 

Turbulence Models 

The turbulence models used in this study axe the two- 
equation £-co model of Wilcox 13,14 , the k-z model of 
Launder and Sharma 15 , the SST (k-w/k-z) model of 
Menter 16 and the algebraic eddy viscosity model formula- 
tion of Baldwin and Lomax 17 . The two-equation k - to 
model is also studied with the compressibility model cor- 
rections, in particular, the length-scale and the separation 
corrections. A detailed description of these models is 
found in the references, and a detailed description of these 
models, including the model corrections is found in refer- 
ences 2 and 4. The main aim of the present effort is to fur- 
ther study the effects of the compressibility corrections in 
the simulation of complex 3-D flows. Currently, research 
on improved compressibility corrections for turbulence 
modeling is being performed to account for the complex 
effects encountered in shock boundary layer interactions. 

Numerical Method 

The cost-effective engineering design of aerospace 
vehicles encountering subsonic, transonic, supersonic and 
hypersonic speeds requires advanced and efficient compu- 
tational fluid dynamics (CFD) technology 8 . Accurate 
aerodynamic prediction of complex full 3-D flow fields 
and the integration of different areas of technology and 
research are presently required to account for the signifi- 
cant nonlinear effects on aerodynamic coefficients, lift, 
drag, and heat load. An improved 3-D Navier-Stokes code 
has been further developed to efficiently validate turbu- 
lence models for high speed flows. The general methodol- 
ogy is found in Bardina 3 , and therefore only a brief 
description of the method is given below. 

The model equations are the 3-D compressible Rey- 
nolds-averaged Navier-Stokes equations Turbulence Mod- 
els in a general curvilinear coordinate system with mass- 
averaged and non-dimensional variables. All flux differ- 
ences are treated implicitly in order to increase stability 
and to be able to use large increments of time or CFL 
numbers. The numerical scheme for the viscous fluxes is 
second-order central difference, while the numerical 
scheme for the in viscid fluxes is a higher-order TVD 
upwind flux-difference splitting. The higher-order TVD 
scheme has the capability to represent first-order upwind, 
second-order upwind, third-order upwind biased, second- 


order Fromm scheme, and other combinations of second- 
order upwind and central differences. 

The efficiency of the method is based on an implicit 
symmetric Gauss-Seidel “method of planes” relaxation 
scheme with alternating directional space marching 
sweeps along one coordinate direction, Newton-Raphson 
inner iteration procedure with an implicit block-tridiago- 
nal diagonally dominant approximate factorization relax- 
ation scheme along the other two directions. This method 
requires less data in central memory and less total transfer 
of data into central memory per iteration than implicit 
upwind schemes using only time-dependent approximate 
factorizations; therefore, the capability of processing 
larger and/or complex data bases and computational grids 
is available. The data is conveniently stored on successive 
planes along the streamwise coordinate, and the system of 
equations is solved twice in each successive plane first 
along the forward direction and afterwards along the 
backward direction. The repeated solution procedure pro- 
vides an effective Newton-Raphson convergence acceler- 
ation. In each plane the solution is obtained by a two level 
diagonally dominant approximate factorization DDADI 
procedure 13 . The space marching alternating directional 
sweeps in the streamwise coordinate are von Neumann 
unconditionally stable for zones of subsonic and stream- 
wise separated and reversed flows as well as supersonic 
flow. As the more restrictive FNS techniques, the present 
space marching method results in improved propagation 
of nonlinear effects to accelerate convergence to steady 
state, generally in about one order of magnitude fewer 
iterations than approximated factorization methods. 

This method combines the best features of data man- 
agement and computational efficiency of space marching 
procedures with the generality and stability of time depen- 
dent Navier-Stokes procedures to solve flows with mixed 
subsonic and supersonic zones, including streamwise sep- 
arated flows. Its robust stability derives from a combina- 
tion of conservative implicit upwind flux difference 
splitting, inner approximation procedure in grid cells 
where changes of eigenvalue sign are present, diagonally 
dominant approximate factorization and relaxation 
scheme, flux limiters of higher-order flux differences, and 
well-posed characteristic-based implicit boundary approx- 
imations. It provides the capability of predicting complex 
flow structures in complex geometries with good accu- 
racy. 

Boundary Conditions 

Mathematically well posed implicit characteristic- 
based boundary procedures were imposed at every bound- 
ary point. The equilibrium turbulent boundary layer was 
prescribed at the inflow boundary points. The inflow pro- 
file matched the experimental displacement thickness, 
5*^=0.0126 m, located at 1.62 m from the leading edge of 
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the flat plate. On the fin and flat plate boundary points, 
constant wall temperature (1^=300 °K) and no slip condi- 
tions were imposed; the turbulent kinetic energy k and its 
dissipation rate e were set equal to zero; (0 was set equal to 
10 times greater than the corresponding theoretical value 
at the first point off the wall. On the symmetry plane, no 
flow through and zero-gradient extrapolation of density, 
pressure, streamwise velocity, and turbulence variables 
were imposed. On the upper free-flow plane and other 
inflow/outflow boundary points, finite difference was 
imposed both along and toward the boundaries. The pro- 
cedure automatically determined whether the fluid was 
flowing locally toward or fromm the boundary, and it 
imposed appropriate conditions accordingly. If the inflow 
was subsonic, no changes in entropy, tangential velocity 
components, enthalpy, and turbulence variables were 
imposed. If the outflow was subsonic, no pressure gradi- 
ent was imposed since only one characteristic-based 
boundary approximation was required in the differences 
toward the boundary. If the outflow was supersonic, the 
solution was naturally extrapolated with the upwind 
scheme with no external boundary approximations. These 
boundary approximations have been proven to be effec- 
tive in previous simulations and free stream has been 
effectively maintained 3 . 

Code Performance 

The numerical simulations in the Cray Y-MP C90 
supercomputer located at NASA Ames Research Center 
performed at a rate of 51 MIPS and 288 MFLOPS. Simu- 
lations studies were done with 101x61x41 and 231x81x81 
grid points. Small differences between the solutions were 
observed in the surface pressure and heat transfer distribu- 
tions, more differences were observed in the flow struc- 
tures. The fine mesh solutions provided the best resolution 
of the turbulence structures, while the less refine mesh 
solutions are considered accurate for engineering pur- 
poses. Most results presented here were obtained with the 
101x61x41 grid and required less than 6 hours of CPU 
time and less than 600 sweeps (or global iterations) to 
achieve convergence to steady state. The performance of 
this diagonal-dominant implicit upwind code shows at 
least one order of magnitude better efficiency than other 
Navier-Stokes codes based on well-known central-differ- 
ence numerical methods. 

Intersecting Shock-Waves/Turbulent Boundary Laver 
Interaction flSWBLI) 

The Ames experiment of Kussoy and Horstman 12 on 
3-D shock-wave boundary-layer interactions was used 
here to test the compressible turbulence models and the 
model corrections. This experiment studies the interac- 
tions of two intersecting hypersonic shock waves with a 
thick turbulent boundary layer. The experimental configu- 
rations reflect several key elements of generic hypersonic 


inlets, thick turbulent boundary-layer approaching two 
vertical fins of varying wedge angles, crossing shock- 
waves, boundary-layer cross-flow vortices, and large pres- 
sure gradients. The test body for this series of experiments 
is shown in fig.l. Two 15° fins mounted on top of a 2.2 m 
long flat plate generated two planar oblique crossing 
shock waves on a thick turbulent boundary layer. The 
free-stream Mach number was Nt = 8.3, the free-stream 
temperature was T^ = 80 °K, the Reynolds number was 
Re^ = 5.3* 10 6 per meter, and the wall temperature was 
fixed at 300 °K. 

The physics of this flow shows a pattern of intersect- 
ing shock-waves above the boundary layer and a complex 
set of cross-flow vortices and structures inside of the tur- 
bulent boundary layer. Previous experimental and compu- 
tational analyses have provided a general description of 
the flow fields generated through the interaction of a sin- 
gle shock-wave and a turbulent boundary layer. Settles 
and Dolling 18 reviewed the early work on this class of tur- 
bulent flow interaction, while Kubota and Stollery 19 
described the main vortical structure developed inside the 
boundary layer and under an oblique shock-wave. The 
interaction of each shock wave with the boundary layer 
generated a cross-flow vortex separation with a “quasi 
conical” shape 1,20 * 21 . Although the “quasi-conicaT struc- 
ture has been used in different studies of turbulence mod- 
els and conical simulations 21,22 , this approximation has 
been disputed previously 1 . The comparison of Bardina et 
al 1 and Knight et al 21 shows that this approximation intro- 
duces large errors and make comparisons of turbulence 
models meaningless. In this particular experiment, the 
influence of the lateral fins in the flow structures and sur- 
face quantities imply the necessity of a full 3-D numerical 
simulation. This flow is further complicated by the inter- 
section of both counter-rotating cross-flow vortices, which 
uplifted the flow, producing large losses of total pressure, 
and generating a very complex flow structure with sec- 
ondary structures developed under the cross-flow vortices 
and on the lateral fins. The fins developed their own 
hypersonic laminar boundary layer with expansion fans 
and lateral separation. In recent numerical investigations, 
Narayaswami et al 23 used the algebraic mixing-length 
model of Baldwin and Lomax, and the modified k-z model 
of Rodi 22 for the turbulent eddy viscosity. Their results 
showed qualitative agreement with experimental data, 
peak surface pressures and heat transfer data are overpre- 
dicted. Gaitonde and Shang 24 have used the Baldwin and 
Lomax turbulence model and Roe’s flux -difference split 
upwind numerical scheme. Their results show agreement 
with surface pressure data, and overprediction of surface 
heat transfer data. 

A selected comparison between experimental data 
and numerical simulation results is described below. In 
general, symbols in the figures shown below represent the 
experimental data points, the solid lines show solutions 
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obtained with two-equation turbulence models, and the dash 
lines show solutions obtained with two-equation turbulence 
models and model corrections (length-scale and rapid com- 
pression corrections). 

Velocity vectors 

A set of velocity vector plots are shown in Figure 2. These 
results were obtained with the jfc-o model including the length- 
scale and rapid compression corrections. The fin boundary lay- 
ers were treated as turbulent below and laminar above the edge 
of the flat-plate boundary layer. Fig. 2a shows the velocity vec- 
tors next to the flat plate surface. It shows the vortex interac- 
tion zone, the flow turning and reflections, and the wake-like 
structure in the downstream zone. Fig. 2b shows the velocity 
vectors in the symmetry plane between the lateral fins. The 
main results show the uplifting of the boundary layer flow due 
to the vortex “collision”, and the secondary uplifting and reat- 
tachment below the main vortices. Fig. 2c, 2d,2e, and 2f show 
velocity vectors in different crossed sections, x/b^ = 3, 6, 9, 
and 12, respectively. They show the formation of two cross- 
flow vortices as main structures, a center bubble under the 
main vortices, comer vortices in the fin/plate junctions, and 
strong flow turnings at the edge of the flat plate boundary layer. 
The strength of the main vortices is model dependent. Figures 
2f, 2g, and 2h compared the standard k-w model with and with- 
out model corrections and modeling the fin boundary layer as 
turbulent/laminar, turbulent, and laminar, respectively. The 
main feature is the increase of vorticity generated by the model 
corrections (Fig. 2f and 2h) and the almost disappearance of 
the main vortices with the standard £-co model (Fig. 2g). 

Pressure contours 

Figures 3a and 3b show the normalized pressure contours 
in two cross-section planes. It shows the well known structure 
of the “quasi-conical” shape, the vortex, the triple point, the 
slip line. In the center zone, a secondary structure on the plate 
surface and a pressure wave between the plate surface and the 
free stream are present. These features show the complexity 
generated by the double fin interaction, beyond the single fin 
case. 

Flow-field yaw angle 

Figures 4a through 4e show the comparison of yaw-angle 
contours between simulation and experiment on different 
cross-sectional planes. The main emphasis is to differentiate 
between the simulations with and without model corrections, 
and to verify if the treatment of the fin boundary layers affect 
the results. The flow turning as represented by the yaw-angle is 
well modeled with the model corrections in Fig. 4b and 4e, 
especially if we consider the few measurements available in 
each plane. On the other hand, the absence of the model correc- 
tions in Fig. 4d shows almost no turning above the main vorti- 
ces as shown in the experimental data. 

Flow-field pitot pressure 

Figures 5a through 45 show the comparison of Pitot pres- 
sure contours between simulation and experiment on different 


cross-sectional planes. The main emphasis is also to differenti- 
ate between the simulations with and without model correc- 
tions The low Pitot pressure zone in the main vortices zone 
shows agreement with the experimental data, especially if we 
consider the few experimental measurement points in each 
plane. The model corrections show larger gradients in Fig. 5b 
than the simulations without the model corrections in Fig. 5d. 
The laminar treatment of the fin boundary layers shown in Fig. 
5e shows large differences with the turbulent/laminar treatment 
shown in Fig. 5b. 

Surface pressure 

Figure 6 shows the comparisons between prediction and 
experimental data of surface pressure distributions on the flat 
plate. Fig. 6a shows the centerline profiles between the two lat- 
eral fins. Figures 6b, 6c, and 6d show transverse surface pres- 
sure distributions at x/8^ = 5.6, 6.92, and 8.31 boundary layer 
thicknesses downstream of the fin leading edges, respectively. 
A wave structure with high and low peaks of static pressure are 
observed. The pressure rises induced by the shock waves gen- 
erated by the fin leading edges, and decreases induced by the 
expansion fans generated by the interior fin corners. The peak 
pressure value is located at the reattachment zone of the sec- 
ondary structure formed under the two principal counter-rotat- 
ing vortices. The first transverse distribution is located in the 
shock-wave “collision” zone, the second one is located near 
the peak surface pressure zone generated by the secondary flow 
re-attachment, and the third one is located near the lower 
expansion pressure zone. The turbulence models shown here 
are the k - co model, the SST model and the Baldwin-Lomax 
model. The k-i o model includes the (1) length-scale, (r) separa- 
tion bubble, and/or (w) rotation corrections. All simulations 
show good agreement with the experimental data, and the peak 
values are well predicted within the experimental uncertainty 
of 10%. Overprediction near the outflow zone is observed 
when the fin boundary layer are computed as (tur) turbulent 
boundary layers. If the fin boundary layers are computed as 
(lam) laminar boundary layers, no overprediction is obtained. 
The best treatment is obtained when the fin boundary layer is 
treated as (tur/lam) turbulent below and laminar above the tur- 
bulent boundary layer on the flat plate. Different from the 
results observed with a single fin 1 , the modeling of the fin 
boundary layers affect the centerline pressure distribution. 
Since the pressure ratio in this 3-D simulations is smaller than 
the ones presented in the 2-D hypersonic database 
experiments 9-11 , the model corrections show only small differ- 
ences in the numerical predictions. 

Surface heat transfer rate 

Figure 7 shows the comparisons between prediction and 
experimental data of surface heat transfer distributions on the 
flat plate. Similar to the pressure distributions shown in the 
previous figures 6, Fig. 7a shows the centerline distributions 
and Figures 7b, 7c, and 7d show the transverse distributions at 
x/S«, = 5.08, 6.4, and 7.78, respectively. A similar wave struc- 
ture with high and low peaks along the centerline is also 
observed. The experimental data in the transverse profiles 
show a flatter distribution in the cross-sectional planes. Recent 
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simulations of Narayanswami et al 23 and Gaitonde et al 24 
showed overprediction of heat transfer rate and “are thought to 
be associated with deficiencies in turbulence modeling.” 
Therefore, the testing of the turbulence model corrections is of 
great interest here. 

All models showed excellent agreement with the experi- 
mental data along the centerline distribution. The small plateau 
at the beginning of the interaction is not shown in this 
101x61x41 grid simulations, however, they are present in the 
231x81x81 simulations not shown here. In the transverse dis- 
tributions, the k-w, SST, and Baldwin-Lomax models show an 
overprediction which seems to be associated with the cross- 
flow reattachment. The model corrections improve the predic- 
tions and show good agreement with the experimental data in 
the first two stations shown in Fig 7b and 7c. This agreement is 
only present when the fin boundary layers are treated as turbu- 
lent ones below the edge of the boundary layer on the flat plate. 
These results support the model corrections and the proper 
treatment of each boundary layer. Once again, best treatment is 
obtained when the fin boundary layer is treated as (tur/lam) tur- 
bulent below and laminar above the turbulent boundary layer 
of the flat plate. In the last transverse station shown in Fig. 7d, 
overprediction in the lateral cross-flow reattachment zone is 
still observed with all models, and it shows the complexity of 
this flow. 

Total pressure 

Figure 8 shows the normalized total pressure distribution 
along the streamwise direction. The numerical results 21 
obtained by Horstman with Rodi’s modified k-t model and 
Knigth with the Baldwin-Lomax model are also included in 
this figure. The results show a significant loss of about 85% in 
total pressure due to the boundary layer interaction with the 
shock-waves forming the two cross-flow vortices as a low total 
pressure outflow jet. The need to eliminate these efficiencies in 
this kind of interaction is a subject of continuous research, 
including boundary layer bleeding and geometry modifica- 
tions. 

Concluding Remarks 

In this section we summarize the research work, give our 
principal results and recommendations, and discuss plans for 
future work.The more promising turbulence model corrections 
for compressible flows were tested. The agreement with the 
experimental data is very good in surface pressure, heat trans- 
fer r rates, yaw angles, and Pitot pressure. The model correc- 
tions give improved heat transfer predictions. Different than 
the single fin simulations. The treatment of the fin boundary 
layer affects the surface plate predictions. The best results are 
obtained with the proper turbulent/laminar boundary layer on 
the fin walls. Accurate and efficient aerodynamic predictions of 
the intersecting SWTBL interaction have been presented. The 
present results show that numerical solutions can be efficiently 
obtained in order to provide a data set for engineering design. 
The flow structures are well captured within a few grid points 
and free of oscillations. The predictions in these zones are 
superior and show detailed primary and secondary turbulence 


structures. The physical understanding of these structures is 
fundamental to improve inlet designs and to improve the com- 
pressibility model corrections. This methodology provides a 
promising computational capability for aerospace vehicles. 
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Fig. 1. Geometry of flat-plate with two 15° lateral 
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Fig. 2c. Velocity vectors on crossflow plane located 
at x/8„ = 3Jfc-co model with model corrections. 
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Fig. 2a. Velocity vectors on top of flat-plate surface, 
fe-co model with model corrections. 


Fig. 2d. Velocity vectors on crossflow plane located 
at x/Seo = 6. )t-(D model with model corrections. 
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Fig. 2b. Velocity vectors on symmetry plane (Z = 0) 
between lateral turbulent/laminar boundary layer fins, k- 
(0 model with model corrections. 


Fig. 2e. Velocity vectors on crossflow plane located 
at x/5„ = 9. k-(0 model with model corrections. 
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Fig. 2f. Velocity vectors on crossflow plane located Fig. 3a. Normalized pressure contours on crossflow 

at xfa. = 12 and turbulent/lam inar fin. k- to model with plane located atJt/6„=6 and turbulent/laminar fin. k- 

model corrections. co model with model corrections. 



Fig. 2g. Velocity vectors on crossflow plane located Fig. 3b. Normalized pressure contours on crossflow 

at x/d m = 12 and turbulent fin. k-(o model with no plane located at jc/8„ = 9 and turbulent/laminar fin. *-co 

model corrections. model with model corrections. 



Fig.2h. Velocity vectors on crossflow plane located Fig. 4a. Yaw angle contours on crossflow plane 

at jc/8„ = 12 and laminar fin. *-co model with model located at x/b^ = 5.6 and turbulent/laminar fin with k-w 

corrections. model and model corrections. 
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Simulation Experiment 
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Fig. 4b. Yaw angle contours on crossflow plane 
located at jc/ 5*, = 6.92 and turbulenl/laminar fin with k- co 
model and model corrections. 



Fig. 4e. Yaw angle contours on crossflow plane 
located at xfb^ = 6.92 and laminar fin with k- o) 
model and model corrections. 


Simulation Experiment 



Fig. 4c. Yaw angle contours on crossflow plane 
located at xfb^ = 8.3 and turbulent/laminar fin with k- 
co model and model corrections. 


Simulation Experiment 
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Fig. 5a. Pitot pressure contours on crossflow 
plane located at x/5„ 5.6 and turbulent/laminar fin 
with k-<a model and model corrections. 


Simulation Experiment 


Simulation Experiment 



Fig. 4d. Yaw angle contours on crossflow plane Fig. 5b. Pitot pressure contours on crossflow 

located at x/8„ = 6.92 and turbulent fin with /fc-to plane located at x/8^ = 6.92 and turbulent/laminar 

model and no model corrections fin with k-V) model and model corrections. 
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Simulation Experiment 



Fig. 5c. Pitot pressure contours on crossflow plane 
located atj/S*, = 8.3 and turbulent/laminar fin with k- 
0) model and model corrections. 



Fig. 6a. Normalized pressure distribution on cen- 
terline plane, z/6*, = 0 and y/8*, = 0. 


Simulation Experiment 
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Fig. 5d. Pitot pressure contours on crossflow plane 
located at jc/ 5„ = 6.92 and turbulent fin with jfc-oi 
model and no model corrections. 


Simulation Experiment 
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Fig. 5e. Pitot pressure contours on crossflow plane 
located at x/8„ = 6.92 and laminar fin with k-i o model 
and model corrections. 



Fig. 6b. Normalized pressure distribution on cross- 
flow plane located at x/8<«, = 5.6 and y/S„ = 0. 
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Fig. 6c. Normalized pressure distribution on cross- 
flow plane located at x/6„ = 6.92 and y/5„ = 0. 
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Fig. 6d. Normalized pressure distribution on cross- 
flow plane located at x/8„ = 8.3 and y/8„ = 0. 
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Fig. 7c. Surface heat transfer distribution on 
crossflow plane located at x/8„ = 6.4 and y/b„ 
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Fig. 7a. Surface heat transfer distribution on center- Fig. 7<L Surface heat transfer distribution on 

line plane, z/8„ = 0 and y/8„ = 0. crossflow plane located at x/5„, = 7.78 and y/8^ = 

0 . 
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Fig. 7b. Surface heat transfer distribution on cross- Fig. 8. Normalized integrated total pressure dis- 

flow plane located at x/8„ = 5.08 and y/8„ =0. tribution. 
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Abstract 

A new computational efficient 3-D compressible Rey- 
nolds-averaged implicit Navier-Stokes method with 
advanced two-equation turbulence models for high speed 
flows is presented. All convective terms are modeled using 
an entropy satisfying higher-order Total Variation Dimin- 
ishing (TVD) scheme based on implicit upwind flux-dif- 
ference split approximations and arithmetic averaging 
procedure of primitive variables. This method combines 
the best features of data management and computational 
efficiency of space marching procedures with the general- 
ity and stability of time dependent Navier-Stokes proce- 
dures to solve flows with mixed supersonic and subsonic 
zones, including stream wise separated flows. Its robust 
stability derives from a combination of conservative 
implicit upwind flux-difference splitting with Roe’s prop- 
erty U to provide accurate shock capturing capability that 
non-conservative schemes do not guarantee, alternating 
symmetric Gauss-Seidel “method of planes” relaxation 
procedure coupled with a three-dimensional two-factor 
diagonal-dominant approximate factorization scheme, 
TVD flux limiters of higher-order flux differences satisfy- 
ing realizability, and well-posed characteristic-based 
implicit boundary-point approximations consistent with 
the local characteristics domain of dependence. The effi- 
ciency of the method is highly increased with Newton- 
Raphson acceleration which allows convergence in essen- 
tially one forward sweep for supersonic flows. 

The method is verified by comparing with experiment 
and other Navier-Stokes methods. Here, results of adia- 
batic and cooled flat plate flows, compression comer flow, 
and 3-D hypersonic shock-wave/turbulent boundary layer 
interaction flows are presented. The robust 3-D method 
achieves a better computational efficiency of at least one 
order of magnitude over the CNS Navier-Stokes code. It 
provides cost-effective aerodynamic predictions in agree- 
ment with experiment, and the capability of predicting 
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complex flow structures in complex geometries with good 
accuracy. 

Introduction 

A new improved 3-D Navier-Stokes code 1,2 has been 
developed to efficiently validate turbulence models for 
high speed flows. The cost-effective engineering design of 
aerospace vehicles encountering subsonic, transonic, 
supersonic and hypersonic speeds requires advanced and 
efficient computational fluid dynamics (CFD) 
technology 3 . Accurate aerodynamic prediction of complex 
full 3-D flow fields and the integration of different areas of 
technology and research are presently required to account 
for the significant nonlinear effects on aerodynamic coeffi- 
cients, lift, drag, and heat load. 

In this last decade, a substantial research effort has 
been conducted on the development of upwind methods 
for the numerical treatment of the convective teims of the 
Navier-Stokes equations 4 ' 18. Three generations of numeri- 
cal methods can be differentiated. The first generation is 
represented by McCormack’s predictor-corrector method 4 
for aerodynamic predictions. The second generation is 
represented by implicit central-difference approximation 
schemes 5 * 6 with dimensionally split approximated factor- 
izations, explicit and implicit artificial viscosity terms, and 
explicit boundary approximations. Although these meth- 
ods have been shown not to be von Neuman uncondition- 
ally stable 14 * 16,17 in the absence of smoothing procedures, 
they have been proven effective in transonic flow-field 
predictions in respective time and space marching imple- 
mentations. In highly nonlinear and rapidly varying flow 
fields, performance has not been as satisfactory. The pre- 
dictions are largely governed by the sensitive upstream 
movements of shocks with increasing Mach number, and 
the computational accuracy and numerical stability 
become of major importance 8 ' 17 . Intimately related to the 
shocks are precursor expansions which are equally diffi- 
cult and important to treat accurately. A final computa- 
tional detail available for improvement is the ability to 
accurately compute flow in the vicinity of contact discon- 
tinuities which appear in streamwise and cross-flow sepa- 
rations, for example. 

The third generation of numerical methods is repre- 
sented by upwind split-difference approximation schemes 
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with diagonaUy-dominant relaxation procedures 16 . Differ- 
ent procedures are available to formulate these methods. 
Major differences are found in flux-vector splitting 
methods 8 ' 11 and flux-difference splitting methods 2,141 . 
The laiger dissipation errors in both inviscid and viscous 
calculations associated with flux-vector splitting methods 
have led to further implementations of flux-difference 
methods 1216 . No major differences are observed between 
solutions of finite- volume and finite-difference procedures 
once convergence is achieved, although there are differ- 
ences in the formulation and/or the location of the cell 
areas of the Jacobian matrix of the generalized coordinate 
transformation 10 - 1 5 - 1 6 . 

Solutions based on first-order difference approxima- 
tions contain too much numerical dissipation. Higher- 
order methods are formulated based on interpolation and/ 
or extrapolation procedures of fluxes or dependent vari- 
ables (MUSCL approach) 11 . Total variation diminishing 
schemes (TVD) are upwind or symmetric high-resolution 
schemes within smooth regions with first-order accuracy 
at local extrema 10 . The aim of adding the least amount of 
numerical dissipation over the physical dissipation to 
avoid spurious oscillations in shock regions is bounded by 
the need to obey the second law of thermodynamics and 
avoid entropy violations. 

The present method represents an advance in the third 
generation of flux-difference splitting methods. It is based 
on the mathematical well-posed discrete-approximation of 
the Jacobian matrix of the flux-difference vectors. 

Governing Equations 

The governing equations of motion are the 3-D com- 
pressible Reynolds-averaged Navier-Stokes and turbu- 
lence model equations for the conservation laws of mass, 
momentum, and energy of a perfect gas with mass-aver- 
aged and non-dimensional variables. The conservation 
laws are written in general curvilinear coordinates £,{•*,) 
and compressed vector notation as 



where repeated sub-indices in any term imply summation 
over the index range (/=1>2,3); sub-index j following a 
comma imply partial differentiation with respect to the 
respective Cartesian spatial coordinate direction x } ; J is the 
Jacobian of the spatial-coordinate transformation 

dx { dx 2 dx 3 

J = ( 2 ) 
and e ijn is the third-order isotropic alternating tensor. 

The conservative dependent-variable vector U and 
flux vectors Fj are written in terms of non-dimensional 
mass-averaged variables. 


-> 

u = 


[- 1 


~ - 


- 

p 


PUj 


0 

P«1 


P U^J + P^J 


T i; 

P«2 

■fj= 

pu 2 Uj + p8 2j 

+ 

X 2; 

p« 3 


pu 2 uj + pb 3j 


X 3, 

e 


(e + p)U: 



p, p and T are 

the fluid density, pressure and 


( 3 ) 


tem- 
perature, respectively; «,• are the Cartesian velocity com- 
ponents; p = pT is the non-dimensional equation of state; 
e = pl(y-l) + OJpuM + pk is the total specific energy; k is 
the turbulent kinetic energy; x# is the stress tensor; q ; is 
the heat flux vector, y is the ratio of specific heats at con- 
stant pressure and volume; and 5^ is the second-order 
isotropic Kronecker delta. The first term of the flux vec- 
tors Fj in Eq. 3 represents the inviscid hyperbolic conser- 
vation law terms, and the second term represents the 
laminar and turbulent viscous stresses and heat fluxes. 


The viscous stress tensor x,y and the heat-flux vector 
q t include both molecular and (Reynolds averaged) turbu- 
lent contributions. These variables are modeled using the 
Boussinesq approximation in terms of an eddy viscosity 
as 

T iy - = 2 (ji + p T ) (S,, - 5 l7 S„„/3) - 2p*8/3 (4) 



where p and pj are the molecular and turbulent (eddy) 
viscosities, Pr = CpP/K and Pr T = 0.9 are the molecular 
and turbulent Prandtl numbers (assuming air), Pr k is a tur- 
bulence model parameter, = («, v + u yV )/2 is the strain rate 
tensor. Both last terms of Eq. 4 and Eq. 5 require the tur- 
bulent kinetic energy k, they are neglected in zero-equa- 
tion turbulence models (for example, Baldwin-Lomax and 
Cebeci-Smith models), and usually included in two-equa- 
tion turbulence models (for example, £-to and k-E models). 

Two-Equation Eddv Viscosity Models 

The turbulent eddy viscosity in two-equation models 
is here expressed in terms of the turbulent kinetic energy, 
k, and either the dissipation rate, e, or the specific dissipa- 
tion rate, to, depending on the model. This expression is 

Pr = <yiiP* 2/e = e = ( 0 * (6) 

where Cp is a modeling constant, and is a damping 
function depending on the specific model used. 

The two transport equations are expressed in terms of 
the generic variables Ar-O (where i^=e or 0=to depending 
on the model) by the following formulas given below 
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The source terms are defined as 


S k = P t -p(ok-L t (8) 

S* = (C«, P t -/ 2 C„ poik + C« p^)(i3//t) +L # (9) 

and the production of turbulent kinetic energy is 

p t = 2Pr - S^5) - 2Sjfik/3 (10) 

The model parameters for Wilcox k-ut model and 
Launder-Sharma k-E model are shown in the following 
Tables 1 and 2. 


Table 1: Wilcox k-w Model, t) ■ co = zjk 
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Table 2: Launder-Sharma k-t Model, D = e = out 


C el = 1.45 

fv. = exp(-3 .4/( 1 +R 7 2 /50)) 

C c2 = 1.92 

f 2 = \-03exp(-R T 2 ) 

o*= 1 

L k = 2p (ch/k/dxj 2 

o e = 1.3 

L t = 2(|ip r /p) i&u/dXndxJ 2 

C„ = 0.09 

R t = pkVpE 


At present, different turbulence models have been 
implemented in the numerical code for validation studies. 
These models include two-equation turbulence models 19 
(*-©, transitional k-w, BSL, SST, k-E) and the algebraic 
Baldwin-Lomax eddy viscosity model 20 . Compressibility 
corrections for the two-equation turbulence models have 
also been incorporated into the code 1,19 . 

Numerical Method 

The present method is a model of the three-dimen- 
sional, implicit, flux-difference split upwind Navier- 
Stokes equations expressed in Favr6-averaged variables 
and generalized curvilinear coordinates. It provides the 
capability of predicting complex flow structures in com- 
plex geometries with good accuracy. 

This method combines the best features of data man- 
agement and computational efficiency of space marching 
procedures with the generality and stability of time depen- 
dent Navier-Stokes procedures to solve flows with mixed 
subsonic and supersonic zones, including streamwise sep- 
arated flows. Its robust stability derives from a combina- 
tion of conservative implicit upwind flux difference 


splitting, inner approximation procedure in grid cells 
where changes of eigenvalue sign are present, three- 
dimensional diagonally dominant approximate factoriza- 
tion and relaxation scheme, flux limiters of higher-order 
flux differences, and well-posed characteristic -based 
implicit boundary approximations. 

The efficiency of the method is based on an implicit 
symmetric Gauss-Seidel “method of planes” relaxation 
scheme with alternating directional space marching 
sweeps along one coordinate direction, Newton-Raphson 
inner iteration procedure, together with an implicit block- 
tridiagonal diagonally dominant approximate factorization 
relaxation scheme along the other two directions. This 
method requires less data in central memory and less total 
transfer of data into central memory per iteration than 
implicit upwind schemes using only time-dependent 
approximate factorizations; therefore, the capability of 
processing larger and/or complex data bases and computa- 
tional grids is available. The data is conveniently stored on 
successive planes along the streamwise coordinate, and 
the system of equations is solved twice in each successive 
plane along the forward direction, firstly, and along the 
backward direction, afterwards. The repeated solution pro- 
cedure provides an effective Newton-Raphson conver- 
gence acceleration. In each plane the solution is obtained 
by a two level diagonally dominant approximate factoriza- 
tion DDADI procedure 16 . The space marching alternating 
directional sweeps in the streamwise coordinate are von 
Neumann unconditionally stable for zones of subsonic and 
streamwise separated and reversed flows as well as super- 
sonic flow. Much as the more restrictive PNS techniques, 
the present space marching method results in improved 
propagation of nonlinear effects to accelerate convergence 
to steady state, generally in about one order of magnitude 
less iterations than two level linearized implicit methods. 
The method requires only a few 2-D cross-flow data 
planes in core at any marching step and thus with fast data 
transfer, such as is available on supercomputers, can effi- 
ciently treat complex problems requiring very large num- 
bers of mesh points. 

Elux-Difference Jacobian Matrix 4 

All flux differences are treated implicitly in order to 
increase stability and to be able to use large increments of 
time or CFL numbers. The differentiation of each flux-dif- 
ference term of the transport equations the inviscid along 
each curvilinear coordinate is defined as a product of a dis- 
crete Jacobian matrix A and the conservative dependent- 
variable vector 8t/, 

8F = A 8U (10) 

For simplicity, the vector notation and the sub-index j are 
omitted. The Jacobian matrix A is decomposed using a 
similarity transformation as 
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A = (MNT)A(MNT) 1 (11) 

where A is a diagonal matrix whose coefficients are the eigen- 
values of the Jacobian matrix, and ( MNT) is the column-eigen- 
vector matrix of A . The M matrix represents the transformation 
between conservative and primitive variable differences. The 
matrix N represents the transformation between Cartesian and 
curvilinear coordinates. The matrix T represents the transfor- 
mation between primitive and characteristic variable differ- 
ences. 

A major difference between most numerical methods is 
the definition of the coefficients of the Jacobian matrix A. All 
well-posed methods converge to the exact Jacobian dF/U in the 
continuous space, however, they have truncation-error differ- 
ences in the discrete space of finite-differences or finite-vol- 
ume. In the discrete space, there are infinite number of matrices 
that obey the identity of Eq. 10. The general solution is 
obtained from 

= A p + cA h (12) 

the sum of a particular solution that obeys Eq. 10 and an homo- 
geneous solution which may be multiplied by any factor c. A 
simple observation is one particular solution given by Roe’s 
method 12 , another particular solution is the one shown below 
in the present method, and the difference between both of them 
is an homogeneous solution. The main difference between all 
these methods is not the identity of Eq. 10, but the differences 
in each coefficient of the Jacobian matrix A, furthermore, they 
represent different flux changes due to particular dependent 
variable differences. Mathematically, the best approximation is 
the one provided by the discrete representation of the partial 
derivative dF/dU , which is close to the one described below. 

In the present method, the arithmetic averaging is used to 
define the matrix A . The conservation property of all flux-dif- 
ference split methods requires the discrete identity 

8 (aft) s=55a + a5ft (12) 

between any two quantities a and ft. The overbar represents a 
monotonous averaging procedure between two extreme quanti- 
ties. In the present method, the overbar represents arithmetic 
averaging (in Roe’s method, p represents geometric averag- 
ing). 

1 0 0 0 ol 1 

1 0 0 0 -u l 

M _ u 2 0 1 0 0 , a/- 1 = -u 2 

w 3 0 0 1 0 — w 3 

S u { u 2 u 3 1 UjUj-S -u l -u 2 -u 3 1 


0 0 0 0 

10 0 0 
0 10 0 (13) 

0 0 10 


"1 u , u , 
where S = Hr 2 

defined as ^ 


and the primitive variable differences are 


M 1 81/ = (5p, p8«,, p8« 2 , p8u 3 ,5(p/(Y- 1))) (14) 


The coordinate transformation matrices for the flux differences 
in the £, coordinate direction are 


10000 10000 

0 *1.6, *i.6, 5 *u 0 0 x' 2 ,6, x ' x(i2 0 

N = 0 *2.6, *2.6, ?i.2 0 ,N l = o *' 2t 6, *'3.6, 0 

0*3.6, *3.6, ^' 1 - 3 0 0 5’u 5*1,2 5*1.3 0 

00001 iP 0 0 0 1_ 


where each vector is normalized as 


£’l.i = 5l.i- / 


*'., 6 = 

The x symbol is used to denote the inverse matrix coefficients 
of the non-orthogonal matrices, i.e., 

riytJY 


■*1 , T ” (*l f 2'~*l,3'(**1.2 # *1,3')) | J 

*^1,2' = (*1,2' ~*1. V (**1, 2'*1,3') ) 1 ” (*1,2'*1, 3') J 

*1.2' = (*2.34 1 .3-^. J 4l.2)(iyM 


which are all equivalent expressions. 

T is a transformation matrix between primitive and characteris- 
tic variables r 

-p 0 0 p/2 p/2 

0 10 0 0 

0 0 1 0 0 


o oof(l-f)-f(l + f) 


0 00 


2 (y- 1) 2(Y-1) 


-l/poo o (y - i)/yp 
oioo o 
ooio o 

0 0 0 Xi ^z i l fl + il 

pc y p V cJ 

o oo ^ ^n fi-i) 

pc yp v cJ 


and the characteristic variable differences (19) are 
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(MNT) • 5U = 


P yp 


\ b£h | f j + <*)§£ 

c v c yp 
^ i, ' 8b / | f j 

C V C ' vn 


The diagonal coefficients of the eigenvalue matrix A 
are obtained from the solution of the characteristic equa- 
tion of A. For the flux-differences in the Ijj coordinate 
direction are 

*■1 = 5l J»j 

x 2 = = E,ijpUj/p ( 20 ) 

^-4 = Xj ~ 0.5 (X,| + A^) + c 

where 

f = JfaJJTfehJ m 


The split of the flux-difference terms follows the sign 
of their respective eigenvalues as follows 

6/ = Aj-ajF* + A j+05 F- (22) 

AF + =A + AU (23) 

A F~ = A AU (24) 

Unconditionally stable implicit methods are constructed 

by forward upwind flux differences with positive eigenval- 
ues and backward flux differences with negative eigenval- 
ues. Splitting of the convective terms of the two-equation 
turbulence models is also done according to the sign of the 
first eigenvalue shown in Eq. 20. 


Spatial higher-order TVD flux-differences in the right- 
hand-side of the convective terms of the conservation law 
equations are defined by using a general “ minmod ” limiter. 
It is based on tests and generalization of the limiters pro- 
posed by Yee and Harten 10 , and Chakravarthy and 
Osher 15 . 

The flux-limited finite-differences of the transport 
equations for the Navier-Stokes equations, Eq. 10, and for 


the two-equation turbulence model equations, Eq. 7, are 
expressed as 

Aj-o. 5 ^ = (1 -0.5 (qj-o .5 - (q/fyo. 5 ) + ) A + A ^ 5 U (25a) 
V.5^ = (1 -°- 5 toj+0.5 * (q/Oj-0.5 D A--A J+a 5 t/ (25b) 

This scheme is first-order if q=0, second-order central dif- 
ference if q=l, and second-order upwind difference if q=r. 
The limiter functions are defined as a linear blending func- 
tion of upwind differences (<{>=- 1) and central differences 
(<(>= 1) as 

q ± = 0.5[(l+^)'P ± + (l-(J.)q ± i ± ] (26) 

with r+ j+0.5 = ">ax{ 0. A^s^/Aj^j/^ } (27) 

r j-0.5 = rnax ( 0, Aj +0 5 F" / Aj^) 5 F ~ } (28) 

where 4* and <; are the flux-limiter functions of the central 
and upwind differences, respectively. For a second-order 
upwind scheme 

rV 1 = max{minmod{ m+a)!*}, minmod^A+b ) } (29) 

with 0<(a,b)<l. This limiter becomes Yeee’s symmetric 
limiter when a=b=0, Roe’s Superbee when a=b=l, and 
Osher-Chakravarthy limiter when a=0 and 0<b<l. To 
avoid entropy violations the factor b is limited to 0<b<0.5, 
otherwise, violations such as oveiprediction of Mach num- 
bers in shock waves may be predicted (Superbee). The end 
result of the flux limiters on the second-order TVD 
upwind scheme is a first-order flux-difference with a factor 
bounded between 0.5 and 1.5. 

Other limiters implemented into the code are the con- 
tinuous van Leer function 11 

= (r t +lr ± l)(l+lr ± l) ( 30 ) 

and Yee’s symmetric limiters 10 

r t s t = minmod{U + j-} (31) 

1*5* = minmod ( 1 j + } +minmod ( 1 } - 1 (32) 

= minmod{ 1 ,r + j - ,0.5(r + +r - ) ) (33) 

In all these limiters, the “ minmod * ’ function is defined as 
the minimum magnitude between positive or negative 
arguments 

minmod (a,b) = sign(aymox[0,min{\a\ t b'sign(a))] (34) 
and returns a value of zero if the product ab < 0 . 
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Relaxation two-factor DDADI 

The numerical method is based on an implicit “ method 
of plane s” symmetric Gauss-Seidell relaxation scheme. 
The data is conveniently stored on successive planes along 
the streamwise coordinate, and the system of equations is 
solved in each successive plane along the forward direc- 
tion, first, and along the backward direction, afterwards. In 
unseparated supersonic and hypersonic flows, the system 
of equations may be solved along the forward direction 
only. In each plane, the solution is obtained by using a 
two-level pseudo-time-dependent relaxation procedure 
based on a diagonally dominant approximate factorization 
DDADI. The space marching method results in improved 
propagation of nonlinear effects to accelerate convergence 
to steady state, much as do the more restrictive PNS tech- 
niques. 

The diagonal dominant approximate factorization of the 
left-hand-side of the transport equations including the 
implicit viscous-diffusion terms leads to the following 
two-factor block tridiagonal equation sequence for the ^ 
plane relaxation method 

[ -04 52 ) + , D, (At'j'YW* = -RHS nj ' +l (35) 

[-(/4 y ) + ,D,(/iy)']-8(/ = -D'5U* (36) 

The diagonally dominant matrix D involves the first-order 
split Jacobian matrices and the Jacobian matrices of the 
viscous terms of all coordinate directions 

D = / + "(A^+A^+A^) (37) 

and the solution is updated from time step n to time step 
n+1 

U n+l = U n + SU n (38) 

Observe that the RHS of equation 35 has an exponent 
rt,n+ 1 because some terms in the streamwise direction are 
already updated at time step n+1 due to the plane relax- 
ation procedure. A Newton-Raphson acceleration proce- 
dure is obtained by solving each plane twice or more times 
in each relaxation sequence. This procedure updates the 
nonlinear coefficients of the Jacobian matrix A and pro- 
duces significant improvement in the propagation of the 
nonlinear waves. 

Boundary Conditions 

Mathematically well posed implicit boundary point 
procedures are imposed in every boundary plane. The 


characteristic-based numerical procedure imposes implicit 
boundary conditions, the code automatically determines 
appropriate numerical boundary approximations based on 
the input parameters of the 3-D program and whether the 
flux-difference splitting indicates that the information is 
propagated from the boundary toward the interior domain. 
If the inflow is supersonic, the equilibrium turbulent 
boundary layer is prescribed. If the inflow is subsonic, no 
variation of entropy, enthalpy, and tangential velocity are 
imposed. If the outflow is supersonic, the solution is com- 
puted naturally with the upwind scheme without imposed 
external boundary conditions. If the outflow is subsonic, 
one characteristic-based boundary approximation is 
required and no pressure variation is imposed in the differ- 
ences toward the boundary. At inflow/outflow boundary 
points, finite difference is done both along and toward the 
boundary, the procedure automatically determines whether 
the fluid is locally flowing inside or outside off the bound- 
ary and imposes appropriate conditions accordingly. On 
cooled walls, prescribed wall temperature and no slip con- 
ditions are imposed; the turbulent kinetic energy k and its 
dissipation rate e are set equal to zero; co is set equal to its 
theoretical value at each first point off the wall boundary 
and equal to 10 times that value at each boundary point. 
On symmetry planes, no flow through and zero-gradient 
extrapolation of density, pressure, streamwise velocity, 
and turbulence variables are imposed. These boundary 
approximations have been proven to be effective in previ- 
ous simulations and free stream has been effectively 
maintained 3 . 

Numerical Convergence 

All simulations were done in the Cray C-90 at NASA 
Ames Research center. 

In the following numerical convergence cases, the free 
stream is assumed to be air obeying the perfect gas law 
with effective specific heat ratio 1.4, molecular Prandtl 
number 0.72, turbulent Prandd number 0.9. The laminar 
viscosity is specified by using Sutherland’s law. 

Fig. 1 shows the comparison of the convergence given 
by the energy residuals between the CNS code 21 and the 
present code in the numerical simulation of a free-stream 
flow and an adiabatic flat plate flow at Mach 4 and Rey- 
nolds number of 610 5 . Convergence to machine accuracy 
is achieved within one order of magnitude of less itera- 
tions with the new method. 

Fig. 2 shows the convergence history of the energy 
residual in the present code for three different flows, a 
Mach 4 and Reynolds number 6 10 6 free-stream flow, a 
Mach 5 and Reynolds number 2.2 10 7 adiabatic flat plate 
flow, and a Mach 8 and Reynolds number 5- 10 6 cooled flat 
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plate flow with T w /T aw =0.3. Convergence is achieved 
within 30 sweeps in the first two cases and within 100 
sweeps in the last case. 



Fig. 1. Comparison of energy residual history 
between CNS code and new RANS code. 



Fig. 2. Energy residual history of new RANS code. 


Mach 5 Adiabatic Flat Plate 

The mean flow conditions of this flow on the adiahatir 
flat plate are Mach 5 and Reynolds number of 2.2- 10 7 per 
unit length. Fig. 3a shows the comparison of compressible 
skin friction and Reynolds number based on momentum 
thickness with the experimental correlation of van Driest 
n. The predictions with the k-w and the modified k-w 
models agree with the theoretical correlation. Fig. 3b 
shows the comparison of the velocity profile at Re e =10 4 
with the numerical results of Coakley et al 19 obtained with 


a well proven 2-D finite-volume flux-difference splitting 
Navier-Stokes code. 



Fig. 3a. Skin friction on adiabatic flat plate with free- 
stream Mach number 5 and Reynolds number 2.2- 10 7 . 



Fig. 3b. Comparison of velocity profile on flat plate. 


Mach 8 In viscid Shock Reflection 

The inviscid shock reflection of a Mach 8 flow is used 
here to test the flux limiters described in Eq. 29. Several 
simulations were run to machine convergence of residuals 
down to 10 14 . Fig. 4a and 4b show the Mach number plot- 
ted as a surface plot on the reflection plane of a 30° shock 
wave, the upper surface correspond to the free stream 
Mach 8, the other two lower surfaces correspond to the 
Mach number levels after the incident shock and reflected 
shock, respectively. The main objective is to test the least 
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diffusive flux limiter (Superbee) and to check the amount 
of numerical dissipation needed to avoid entropy viola- 
tions. 



Fig. 4a. Mach number with a=l and b=l in Eq. 29. 



Fig. 4b. Mach number with a=l and b=0.5 in Eq. 29. 

The results with the least diffusive Superbee scheme 
shown in Fig. 4a develops local Mach number peaks 
upstream of the incident shock wave. If the amount of sec- 
ond-order central-difference is increased, as it is shown in 
Fig. 4b, these local peaks vanish and the solution is 
smooth. The test show that b should be no larger than 0.5 
in Eq. 29, otherwise entropy violations will occur. This 
effect is very important in turbulent flows because the 
model variables, k and co or e, also develop strong fluctua- 
tions in the presence of local Mach number peaks. This is 
an effect observed in 3-D simulations of inlets. 

Mach 2.84 24° C nmnression Corner 

The experimental data on the separated turbulent 
boundary layer of Settles et al 22 ’ 23 on a 24° compression 
comer with a free-stream Mach number of 2.84 is com- 


pared with numerical results obtained with the present 
method and the k - to tuibulence model. 



Fig. 5a. Comparison of surface pressure on 24° com- 
pression comer flow with free-stream Mach number 2.84. 



Fig. 5b. Skin friction comparison. 

The agreement between the numerical results and the 
experimental data is very good with small underprediction 
of separation due to the turbulence model. 

Mach 8.3 Intersecting Sho ck- Waves/Turbulent Bound- 
ary Laver Interaction 

Figs. 6 show the comparison of a complex 3-D inter- 
secting shock-waves/turbulent boundary layer interaction 
flow with free stream Mach number 8.3, free stream Rey- 
nolds number 510 6 , and wall temperature 0.3 times the 
adia harir. wall temperature. The experimental data of Kus- 
soy et al 24 belongs to the compressible database developed 
for model validation 22 . The experiment consists of a tur- 
bulent boundary layer on a flat plate with two 15° fins gen- 
erating two oblique shocks with cross-flow separation 
within the boundary layer and a very complex vortex 
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interaction. Figs. 6a and 6b show the surface pressure and 
the heat transfer profiles, respectively, on the flat plate 
along the symmetry plane between the two fins. Fig. 6c 
and 6d show the surface pressure and heat transfer profiles 
on a crossed plane 7 boundary layer thicknesses down- 
stream of the fin edges. The numerical predictions 
obtained with different two-equation turbulence models 
and the algebraic Baldwin-Lomax eddy viscosity model 
are compared with experimental data showing excellent 
agreement with experimental data. 



Fig. 6a. Normalized pressure distribution on center- 
line plane, z/5^ = 0 and y/JL = 0. 



Fig. 6b. Surface heat transfer distribution on center- 
line plane, z/8^ = 0 and y/8^ = 0. 

The centerline profiles show an increase in surface pres- 
sure and heat transfer in the collision zone between the 
cross-flow vortices generated under the oblique shock 
waves. Peak values are obtained in the reattachment zone 
behind the collision zone. Lateral expansion waves 


decrease the magnitude of the profiles, and a wave struc- 
ture develops due to the intersecting shocks. 



Fig. 6c. Normalized pressure distribution on cross- 
section plane located at x/5^ = 6.92 and y/6^ = 0. 



Fig. 6d. Surface heat transfer distribution on cross- 
section plane located at x/5^ = 6.4 and y/8 m =0. 

Concluding Remarks 

Accurate and efficient aerodynamic numerical simu- 
lations have been presented. The results show the robust- 
ness and stability of the 3-D method and code including 
different two-equation turbulence models. Validation is 
shown with predictions of turbulent flows on a Mach 5 
adiabatic flat plate, Mach 3 compression comer, and 3-D 
Mach 8 intersecting shock-waves/turbulent boundary layer 
interaction. The results are in quite good agreement with a 
well validated 2-D code, correlations, and experiment. 

The present results show that numerical solutions of 
turbulent high speed flows can be efficiently obtained in 
order to provide a data set for engineering design. The 
flow structures are generally well captured within a few 
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grid points and are free of oscillations, especially the 
oblique shocks, as well as, the expansion and compression 
waves in the different complex flow zones. This method 
provides a promising computational capability to simulate 
turbulent flows in complex flow fields of aerospace vehi- 
cles. 
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